Potassium humate supplementation improves photosynthesis and agronomic and yield traits of foxtail millet.

Foxtail millet is a highly nutritious crop, which is widely cultivated in arid and semi-arid areas worldwide. Humic acid (HA), as a common plant growth regulator, is used as an organic fertilizer and feed additive in agricultural production. However, the impact of potassium humate KH on the photosynthetic rate and yield of foxtail millet has not yet been studied. We explored the effects of KH application on the morphology, photosynthetic ability, carbon and nitrogen metabolism, and yield of foxtail millet. A field experiment was performed using six concentrations of KH (0, 20, 40, 80, 160, and 320 kg ha-1) supplied foliarly at the booting stage in Zhangza 10 cultivar (a widely grown high-yield variety). The results showed that KH treatment increased growth, chlorophyll content (SPAD), photosynthetic rate (Pn), transpiration rate (Tr), and stomatal conductance (Gs). In addition, soluble protein content, sugar content, and nitrate reductase activity increased in KH-treated plants. With increased KH concentration, the effects became more evident and the peak values of each factor were achieved at 80 kg ha-1. Photosynthetic rate showed significant correlation with SPAD, Tr, Gs, and soluble protein content, but was negatively correlated with intercellular CO2 concentration. Compared to that of the control, the yield of foxtail millet under the T2, T3, T4, and T5 (40, 80, 160, and 320 kg ha-1 of KH) treatments significantly increased by 6.0%, 12.7%, 10.5%, and 8.6%, respectively. Yield exhibited a significant positive correlation with Tr, Pn, and Gs. Overall, KH enhances photosynthetic rate and yield of foxtail millet, therefore it may be conducive to stable millet production. These findings may provide a theoretical basis for the green and efficient production of millet fields.

Potential global distribution of Setaria italica, an important species for dryland agriculture in the context of climate change.

Setaria italica (S. italica, Linnaeus, 1753) is a drought-resistant, barren-tolerant, and widely adapted C-4 crop that plays a vital role in maintaining agricultural and economic stability in arid and barren regions of the world. However, the potential habitat of S. italica under current and future climate scenarios remains to be explored. Predicting the potential global geographic distribution of S. italica and clarifying its ecological requirements can help promote sustainable agriculture, which is crucial for addressing the global food crisis. In this study, we predicted the potential global geographic distribution of S. italica based on 3,154 global distribution records using the Maxent model and ArcGIS software. We assessed the constraints on its potential distribution based on the contribution of environmental factors variables. The predictive accuracy of the Maxent model was evaluated using AUC values, TSS values, and Kappa statistics, respectively. The results showed that the Maxent model had a high prediction accuracy, and the simulation results were also reliable; the total suitable habitats of S. italica is 5.54×107 km2, which mainly included the United States (North America), Brazil (South America), Australia (Oceania), China, India (Asia), and the Russian Federation (Europe). The most suitable habitat of S. italica was 0.52×107 km2, accounting for 9.44% of the total areas, mainly in the United States, India, the Russian Federation, and China. Soil and precipitation (driest monthly precipitation, hottest seasonal precipitation) are the most critical factors limiting the potential distribution of S. italica. Compared with the modern potential distribution, we predict that the four future climate change scenarios will result in varying reductions in the possible geographic ranges of S. italica. Overall, climate change may significantly affect the global distribution of S. italica, altering its worldwide production and trade patterns.

Reactive oxygen species (ROS) modulate nitrogen signaling using temporal transcriptome analysis in foxtail millet.

Reactive oxygen species (ROS) is a chemically reactive chemical substance containing oxygen and a natural by-product of normal oxygen metabolism. Excessive ROS affect the growth process of crops, which will lead to the decrease of yield. Nitrogen, as a critical nutrient element in plants and plays a vital role in plant growth and crop production. Nitrate is the primary nitrogen source available to plants in agricultural soil and various natural environments. However, the molecular mechanism of ROS-nitrate crosstalk is still unclear. In this study, we used the foxtail millet (Setaria italica L.) as the material to figure it out. Here, we show that excessive NaCl inhibits nitrate-promoted plant growth and nitrogen use efficiency (NUE). NaCl induces ROS accumulation in roots, and ROS inhibits nitrate-induced gene expression in a short time. Surprisingly, low concentration ROS slight promotes and high concentration of ROS inhibits foxtail millet growth under long-term H2O2 treatment. These results may open a new perspective for further exploration of ROS-nitrate signaling pathway in plants.

Leaf Senescence Regulation Mechanism Based on Comparative Transcriptome Analysis in Foxtail Millet.

Leaf senescence, a pivotal process in plants, directly influences both crop yield and nutritional quality. Foxtail millet (Setaria italica) is a C4 model crop renowned for its exceptional nutritional value and stress tolerance characteristics. However, there is a lack of research on the identification of senescence-associated genes (SAGs) and the underlying molecular regulatory mechanisms governing this process. In this study, a dark-induced senescence (DIS) experimental system was applied to investigate the extensive physiological and transcriptomic changes in two foxtail millet varieties with different degrees of leaf senescence. The physiological and biochemical indices revealed that the light senescence (LS) variety exhibited a delayed senescence phenotype, whereas the severe senescence (SS) variety exhibited an accelerated senescence phenotype. The most evident differences in gene expression profiles between these two varieties during DIS included photosynthesis, chlorophyll, and lipid metabolism. Comparative transcriptome analysis further revealed a significant up-regulation of genes related to polysaccharide and calcium ion binding, nitrogen utilization, defense response, and malate metabolism in LS. In contrast, the expression of genes associated with redox homeostasis, carbohydrate metabolism, lipid homeostasis, and hormone signaling was significantly altered in SS. Through WGCNA and RT-qPCR analyses, we identified three SAGs that exhibit potential negative regulation towards dark-induced leaf senescence in foxtail millet. This study establishes the foundation for a further comprehensive examination of the regulatory network governing leaf senescence and provides potential genetic resources for manipulating senescence in foxtail millet.

Genetic diversity of grain yield traits and identification of a grain weight gene SiTGW6 in foxtail millet.

KEY MESSAGE: Agronomic traits were evaluated in 1250 foxtail millet accessions, and a crucial gene SiTGW6 governing grain yield was identified. Elite haplotypes and dCAPS markers developed for SiTGW6 facilitate molecular breeding. A comprehensive evaluation of phenotypic characteristics and genetic diversity in germplasm resources are important for gene discovery and breeding improvements. In this study, we conducted a comprehensive evaluation of 1250 foxtail millet varieties, assessing seven grain yield-related traits and fourteen common agronomic traits over two years. Principal component analysis, correlation analysis, and cluster analysis revealed a strong positive correlation between 1000-grain weight and grain width with grain yield, emphasizing their importance in foxtail millet breeding. Additionally, we found that panicle weight positively correlated with 1000-grain weight but negatively correlated with branch and tiller numbers, indicating selection factors during domestication and breeding. Using this information, we identified 27 germplasm resources suitable for high-yield foxtail millet breeding. Furthermore, through an integration of haplotype variations and phenotype association analysis, we pinpointed a crucial gene, SiTGW6, responsible for governing grain yield in foxtail millet. SiTGW6 encodes an IAA-glucose hydrolase, primarily localized in the cytoplasm and predominantly expressed in flowering panicles. Employing RNAseq analysis, we identified 1439 differentially expressed genes across various SiTGW6 haplotypes. Functional enrichment analysis indicating that SiTGW6 regulates grain yield through the orchestration of auxin and glucan metabolism, as well as plant hormone signaling pathways. Additionally, we have identified elite haplotypes and developed dCAPS markers for SiTGW6, providing valuable technical tools to facilitate molecular breeding efforts in foxtail millet.

Physiological Investigation and Transcriptome Analysis Reveals the Mechanisms of Setaria italica's Yield Formation under Heat Stress.

Setaria italica is an important crop in China that plays a vital role in the Chinese dietary structure. In the last several decades, high temperature has become the most severe climate issue in the world, which causes great harm to the yield and quality formation of millet. In this study, two main cultivated varieties (ZG2 and AI88) were used to explore the photosynthesis and yield index of the whole plant under heat stress. Results implied that photosynthesis was not inhibited during the heat stress, and that the imbalance in sugar transport between different tissues may be the main factor that affects yield formation. In addition, the expression levels of seven SiSUT and twenty-four SiSWEET members were explored. Sugar transporters were heavily affected during the heat stress. The expression of SiSWEET13a was inhibited by heat stress in the stems, which may play a vital role in sugar transport between different tissues. These results provide new insights into the yield formation of crops under heat stress, which will provide guidance to crop breeding and cultivation.

Genome-wide identification and expression profiling of the ABI5 gene family in foxtail millet (Setaria italica).

BACKGROUND: ABA Insensitive 5 (ABI5) is a basic leucine zipper transcription factor that crucially influences plant growth, development, and stress response. However, there is minimal research on the ABI5 family in foxtail millet. RESULTS: In this study, 16 ABI5 genes were identified in foxtail millet, and their sequence composition, gene structures, cis-acting elements, chromosome positions, and gene replication events were analyzed. To more thoroughly evaluate the developmental mechanisms of the SiABI5 family during evolution, we selected three dicotyledons (S. lycopersicum, A. thaliana, F. tataricum) and three (Z. mays, O. sativa, S. bicolor) specific representative monocotyledons associated with foxtail millet for comparative homology mapping. The results showed that foxtail millet ABI5 genes had the best homology with maize. A promoter sequence analysis showed that the SiABI5s contain numerous cis-acting elements related to hormone and stress responses, indicating that the regulation of SiABI5 expression was complex. The expression responses of 16 genes in different tissues, seed germination, and ear development were analyzed. A total of six representative genes were targeted from five subfamilies to characterize their gene expression responses to four different abiotic stresses. Overexpression of SiABI5.12 confers tolerance to osmotic stress in transgenic Arabidopsis thaliana, which demonstrated the function of SiABI5 responded to abiotic stress. CONCLUSIONS: In summary, our research results comprehensively characterized the SiABI5 family and can provide a valuable reference for demonstrating the role of SiABI5s in regulating abiotic stress responses in foxtail millet.

Genome-wide identification of AAAP gene family and expression analysis in response to saline-alkali stress in foxtail millet (Setaria italica L.).

Amino acid/auxin permease (AAAP) genes encode a large family of protein transporters that play important roles in various aspects of plant growth and development. Here, we performed genome-wide identification of members in the foxtail millet (Setaria italica L.) AAAP family (SiAAAP) and their saline-alkali stress-induced expression patterns, resulting in the identification of 65 SiAAAP genes, which could be divided into eight subfamilies. Except for SiAAAP65, the remaining 64 genes were located on nine chromosomes of foxtail millet. Gene structure and conserved motif analyses indicated that the members in the same subfamily are highly conserved. Gene duplication event analysis suggested that tandem duplication may be the main factor driving the expansion of this gene family, and Ka/Ks analysis indicated that all the duplicated genes have undergone purifying selection. Transcriptome analysis showed differential expression of SiAAAPs in roots, stems, leaves, and tassel inflorescence. Analysis of cis-acting elements in the promoter indicated that SiAAAPs contain stress-responsive cis-acting elements. Under saline-alkali stress, qRT-PCR analysis showed that SiAAP3, SiLHT2, and SiAAP16 were differentially expressed between salt-alkali tolerant millet variety JK3 and salt-alkali sensitive millet variety B175. These results suggest that these genes may be involved in or regulate the response to saline-alkali stress, providing a theoretical basis for further studying the function of SiAAAPs.

The bHLH transcription factor OsPRI1 activates the Setaria viridis PEPC1 promoter in rice.

Photosynthetic efficiency is reduced by the dual role of Rubisco, which acts either as a carboxylase or as an oxygenase, the latter leading to photorespiration. C4 photosynthesis evolved as a carbon-concentrating mechanism to reduce photorespiration. To engineer C4 into a C3 plant, it is essential to understand how C4 genes, such as phosphoenolpyruvate carboxylase (PEPC1), are regulated to be expressed at high levels and in a cell-specific manner. Yeast one-hybrid screening was used to show that OsPRI1, a rice bHLH transcription factor involved in iron homeostasis, binds to the Setaria viridis PEPC1 promoter. This promoter drives mesophyll-specific gene expression in rice. The role of OsPRI1 in planta was characterized using a rice line harbouring SvPEPC1pro ::GUS. We show that OsPRI1 activates the S. viridis PEPC1 promoter by binding to an N-box in the proximal promoter, and that GUS activity is highly reduced in SvPEPC1pro ::GUS lines when OsPRI1 is mutated. Cross-species comparisons showed that the SvPRI1 homolog binds to the SvPEPC1 promoter but the maize ZmPRI1 does not bind to the ZmPEPC1 promoter. Our results suggest that elements of the iron homeostasis pathway were co-opted to regulate PEPC1 gene expression during the evolution of some but not all C4 species.

Effects of Setaria viridis on heavy metal enrichment tolerance and bacterial community establishment in high-sulfur coal gangue.

Plant enrichment and tolerance to heavy metals are crucial for the phytoremediation of coal gangue mountain. However, understanding of how plants mobilize and tolerate heavy metals in coal gangue is limited. This study conducted potted experiments using Setaria viridis as a pioneer remediation plant to evaluate its tolerance to coal gangue, its mobilization and enrichment of metals, and its impact on the soil environment. Results showed that the addition of 40% gangue enhanced plant metal and oxidative stress resistance, thereby promoting plant growth. However, over 80% of the gangue inhibited the chlorophyll content, photoelectron conduction rate, and biomass of S. viridis, leading to cellular peroxidative stress. An analysis of metal resistance showed that endogenous S in coal gangue promoted the accumulation of glutathione, plant metal chelators, and non-protein thiols, thereby enhancing its resistance to metal stress. Setaria viridis cultivation affected soil properties by decreasing nitrogen, phosphorus, conductivity, and urease and increasing sucrase and acid phosphatase in the rhizosphere soil. In addition, S. viridis planting increased V, Cr, Ni, As, and Zn in the exchangeable and carbonate-bound states within the gangue, effectively enriching Cd, Cr, Fe, S, U, Cu, and V. The increased mobility of Cd and Pb was correlated with a higher abundance of Proteobacteria and Acidobacteria. Heavy metals, such as As, Fe, V, Mn, Ni, and Cu, along with environmental factors, including total nitrogen, total phosphorus, urease, and acid phosphatase, were the primary regulatory factors for Sphingomonas, Gemmatimonas, and Bryobacter. In summary, S. viridis adapted to gangue stress by modulating antioxidant and elemental enrichment systems and regulating the release and uptake of heavy metals through enhanced bacterial abundance and the recruitment of gangue-tolerant bacteria. These findings highlight the potential of S. viridis for plant enrichment in coal gangue areas and will aid the restoration and remediation of these environments.

Efficient identification of QTL for agronomic traits in foxtail millet (Setaria italica) using RTM- and MLM-GWAS.

KEY MESSAGE: Eleven QTLs for agronomic traits were identified by RTM- and MLM-GWAS, putative candidate genes were predicted and two markers for grain weight were developed and validated. Foxtail millet (Setaria italica), the second most cultivated millet crop after pearl millet, is an important grain crop in arid regions. Seven agronomic traits of 408 diverse foxtail millet accessions from 15 provinces in China were evaluated in three environments. They were clustered into two divergent groups based on genotypic data using ADMIXTURE, which was highly consistent with their geographical distribution. Two models for genome-wide association studies (GWAS), namely restricted two-stage multi-locus multi-allele (RTM)-GWAS and mixed linear model (MLM)-GWAS, were used to dissect the genetic architecture of the agronomic traits based on 13,723 SNPs. Eleven quantitative trait loci (QTLs) for seven traits were identified using two models (RTM- and MLM-GWAS). Among them, five were considered stable QTLs that were identified in at least two environments using MLM-GWAS. One putative candidate gene (SETIT_006045mg, Chr4: 744,701-746,852) that can enhance grain weight per panicle was identified based on homologous gene comparison and gene expression analysis and was validated by haplotype analysis of 330 accessions with high-depth (10×) resequencing data (unpublished). In addition, homologous gene comparison and haplotype analysis identified one putative foxtail millet ortholog (SETIT_032906mg, Chr2: 5,020,600-5,029,771) with rice affecting the target traits. Two markers (cGWP6045 and kTGW2906) were developed and validated and can be used for marker-assisted selection of foxtail millet with high grain weight. The results provide a fundamental resource for foxtail millet genetic research and breeding and demonstrate the power of integrating RTM- and MLM-GWAS approaches as a complementary strategy for investigating complex traits in foxtail millet.

Transcriptome analysis and identification of the low potassium stress-responsive gene SiSnRK2.6 in foxtail millet (Setaria italica L.).

KEY MESSAGE: The transcriptome is beneficial for dissecting the mechanism of millet in response to low potassium stress and SiSnRK2.6 was identified as a potential target for improving low potassium stress tolerance. Foxtail millet (Setaria italica L.), which originated in China, has high nutrient utilization character. Nevertheless, the molecular mechanism of its tolerance to low potassium stress is largely unclear. In this research, the low potassium tolerant variety "Yugu28" was screened out by low potassium stress treatment, and the transcriptome of "Yugu28" under low potassium stress was comprehensively analyzed. A total of 4254 differentially expressed genes (DEGs) were identified, including 1618 up-regulated and 2636 down-regulated genes, respectively. In addition, there were 302 transcription factor (TF) genes in the DEGs and MYB TFs accounted for the highest proportion, which was 14.9%. After functional analysis of all DEGs, a total of 7 genes involved in potassium transport and potassium ion channels and 50 genes corresponding to hormones were screened. The expression levels of randomly selected 17 DEGs were verified by qRT-PCR and the results coincided well with the RNA-seq analysis, indicating the reliability of our transcriptome data. Moreover, one of the ABA signaling pathway genes, SiSnRK2.6, was identified and selected for further functional verification. Compared with the wild type, transgenic rice with ecotopic expression of SiSnRK2.6 showed remarkably increased root length and root number, indicating that overexpression of SiSnRK2.6 can enhance the resistance of transgenic plants to low potassium stress.

Comparative transcriptome analysis and genetic dissection of vegetative branching traits in foxtail millet (Setaria italica).

KEY MESSAGE: Two major genetic loci, qTN5.1 and qAB9.1, were identified and finely mapped to the 255 Kb region with one potential candidate gene for tiller number and the 521 Kb region with eight candidate genes for axillary branch number, respectively. Vegetative branching including tillering and axillary branching are vital traits affecting both the plant architecture and the biomass in cereal crops. However, the mechanism underlying the formation of vegetative branching in foxtail millet is largely unknown. Here, a foxtail millet cultivar and its bushy wild relative Setaria viridis accession were used to construct segregating populations to identify candidate genes regulating tiller number and axillary branch number. Transcriptome analysis using vegetative branching bud samples of parental accessions was performed, and key differentially expressed genes and pathways regulating vegetative branching were pointed out. Bulk segregant analysis on their F2:3 segregating population was carried out, and a major QTL for tiller number (qTN5.1) and two major QTLs for axillary branch number (qAB2.1 and qAB9.1) were detected. Fine-mapping strategy was further performed on F2:4 segregate population, and Seita.5G356600 encoding ß-glucosidase 11 was identified as the promising candidate gene for qTN5.1, and eight genes, especially Seita.9G125300 and Seita.9G125400 annotated as B-S glucosidase 44, were finally identified as candidate genes for regulating axillary branching. Findings in this study will help to elucidate the genetic basis of the vegetative branching formation of foxtail millet and lay a foundation for breeding foxtail millet varieties with ideal vegetative branching numbers.

GOLDEN2-like1 is sufficient but not necessary for chloroplast biogenesis in mesophyll cells of C4 grasses.

Chloroplasts are the site of photosynthesis. In land plants, chloroplast biogenesis is regulated by a family of transcription factors named GOLDEN2-like (GLK). In C4 grasses, it has been hypothesized that genome duplication events led to the sub-functionalization of GLK paralogs (GLK1 and GLK2) to control chloroplast biogenesis in two distinct cell types: mesophyll and bundle sheath cells. Although previous characterization of golden2 (g2) mutants in maize has demonstrated a role for GLK2 paralogs in regulating chloroplast biogenesis in bundle sheath cells, the function of GLK1 has remained elusive. Here we show that, contrary to expectations, GLK1 is not required for chloroplast biogenesis in mesophyll cells of maize. Comparisons between maize and Setaria viridis, which represent two independent C4 origins within the Poales, further show that the role of GLK paralogs in controlling chloroplast biogenesis in mesophyll and bundle sheath cells differs between species. Despite these differences, complementation analysis revealed that GLK1 and GLK2 genes from maize are both sufficient to restore functional chloroplast development in mesophyll and bundle sheath cells of S. viridis mutants. Collectively our results suggest an evolutionary trajectory in C4 grasses whereby both orthologs retained the ability to induce chloroplast biogenesis but GLK2 adopted a more prominent developmental role, particularly in relation to chloroplast activation in bundle sheath cells.

Histone deacetylase 9 interacts with SiHAT3.1 and SiHDA19 to repress dehydration responses through H3K9 deacetylation in foxtail millet.

Climate change inflicts several stresses on plants, of which dehydration stress severely affects growth and productivity. C4 plants possess better adaptability to dehydration stress; however, the role of epigenetic modifications underlying this trait is unclear. In particular, the molecular links between histone modifiers and their regulation remain elusive. In this study, genome-wide H3K9 acetylation (H3K9ac) enrichment using ChIP-sequencing was performed in two foxtail millet cultivars with contrasting dehydration tolerances (IC403579, cv. IC4-tolerant, and IC480117, cv. IC41-sensitive). It revealed that a histone deacetylase, SiHDA9, was significantly up-regulated in the sensitive cultivar. Further characterization indicated that SiHDA9 interacts with SiHAT3.1 and SiHDA19 to form a repressor complex. SiHDA9 might be recruited through the SiHAT3.1 recognition sequence onto the upstream of dehydration-responsive genes to decrease H3K9 acetylation levels. The silencing of SiHDA9 resulted in the up-regulation of crucial genes, namely, SiRAB18, SiRAP2.4, SiP5CS2, SiRD22, SiPIP1;4, and SiLHCB2.3, which imparted dehydration tolerance in the sensitive cultivar (IC41). Overall, the study provides mechanistic insights into SiHDA9-mediated regulation of dehydration stress response in foxtail millet.

Underlying Mechanisms of the Hedgehog-Like Panicle and Filamentous Leaf Tissue Symptoms Caused by Sclerospora graminicola in Foxtail Millet.

Downy mildew caused by Sclerospora graminicola is a systemic infectious disease affecting foxtail millet production in Africa and Asia. S. graminicola-infected leaves could be decomposed to a state where only the veins remain, resulting in a filamentous leaf tissue symptom. The aim of the present study was to investigate how S. graminicola influences the formation of the filamentous leaf tissue symptoms in hosts at the morphological and molecular levels. We discovered that vegetative hyphae expanded rapidly, with high biomass accumulated at the early stages of S. graminicola infection. In addition, S. graminicola could affect spikelet morphological development at the panicle branch differentiation stage to the pistil and stamen differentiation stage by interfering with hormones and nutrient metabolism in the host, resulting in hedgehog-like panicle symptoms. S. graminicola could acquire high amounts of nutrients from host tissues through secretion of ß-glucosidase, endoglucanase, and pectic enzyme, and destroyed host mesophyll cells by mechanical pressure caused by rapid expansion of hyphae. At the later stages, S. graminicola could rapidly complete sexual reproduction through tryptophan, fatty acid, starch, and sucrose metabolism and subsequently produce numerous oospores. Oospore proliferation and development further damage host leaves via mechanical pressure, resulting in a large number of degraded and extinct mesophyll cells and, subsequently, malformed leaves with only veins left, that is, "filamentous leaf tissue." Our study revealed the S. graminicola expansion characteristics from its asexual to sexual development stages, and the potential mechanisms via which the destructive effects of S. graminicola on hosts occur at different growth stages.

Multi-environment GWAS identifies genomic regions underlying grain nutrient traits in foxtail millet (Setaria italica).

KEY MESSAGE: A total of 104 foxtail millet accessions were evaluated for 11 nutrients in three environments and 67 high-confidence marker-trait associations (MTAs) were identified. Six SNPs showed pleiotropic effect and associated with two or more nutrients, whereas 24 candidate genes were identified for 28 MTAs involving seven traits. Millets are known for their better nutritional profiles compared to major cereals. Foxtail millet (Setaria italica) is rich in nutrients essential to circumvent malnutrition and hidden hunger. However, the genetic determinants underlying this trait remain elusive. In this context, we evaluated 104 diverse foxtail millet accessions in three different environments (E1, E2, and E3) for 11 nutrients and genotyped with 30K SNPs. The genome-wide association study showed 67 high-confidence (Bonferroni-corrected) marker-trait associations (MTAs) for the nutrients except for phosphorus. Six pleiotropic SNPs were also identified, which were associated with two or more nutrients. Around 24 candidate genes (CGs) were identified for 28 MTAs involving seven nutrients. A total of 17 associated SNPs were present within the gene region, and five (5) were mapped in the exon of the CGs. Significant SNPs, desirable alleles and CGs identified in the present study will be useful in breeding programmes for trait improvement.

Comprehensive physiological, transcriptomic, and metabolomic analysis of the key metabolic pathways in millet seedling adaptation to drought stress.

Drought is one of the leading environmental constraints that affect the growth and development of plants and, ultimately, their yield and quality. Foxtail millet (Setaria italica) is a natural stress-resistant plant and an ideal model for studying plant drought resistance. In this study, two varieties of foxtail millet with different levels of drought resistance were used as the experimental material. The soil weighing method was used to simulate drought stress, and the differences in growth, photosynthetic physiology, metabolite metabolism, and gene transcriptional expression under drought stress were compared and analyzed. We aimed to determine the physiological and key metabolic regulation pathways of the drought-tolerant millet in resistance to drought stress. The results showed that drought-tolerant millet exhibited relatively stable growth and photosynthetic parameters under drought stress while maintaining a relatively stable level of photosynthetic pigments. The metabolomic, transcriptomic, and gene co-expression network analysis confirmed that the key to adaptation to drought by millet was to enhance lignin metabolism, promote the metabolism of fatty acids to be transformed into cutin and wax, and improve ascorbic acid circulation. These findings provided new insights into the metabolic regulatory network of millet adaptation to drought stress.